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Abstract 

The mode-dependent transmission of relativistic ballistic massless Dirac fermion through a 
graphene based double barrier structure is being investigated for various barrier parameters. We 
compare our results with already published work and point out the relevance of these findings to a 
systematic study of the transport properties in double barrier structures. An interesting situation 
arises when we set the potential in the leads to zero, then our 2D problem reduces effectively to a 
ID massive Dirac equation with an effective mass proportional to the quantized wave number along 
the transverse direction. Furthermore we have shown that the minimal conductivity and maximal 
Fano factor remain insensitive to the ratio between the two potentials ( = a 
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1 Introduction 



The physics of graphene material was celebrated by attributing the 2010 physics Nobel prize to two 
physicists, Novoselov and Geim, for their pioneering work on graphene. This material was described 
by some physicists as "heavenly" due to its marvelous physical and transport properties [T]. The story 
started about seven years ago when the possibility to isolate and investigate graphene [H O [3]> i.e. 
individual layers of graphite only one-atom-thick, have been demonstrated. Further, experiments were 
reported showing that charge carriers in graphene behave as two-dimensional relativistic particles with 
zero effective mass. In fact, one of the most interesting aspects of the graphene problem is that its 
low-energy excitations are massless, chiral, Dirac fermions. The quasiparticle excitations around the 
Dirac point obey a linear type of energy dispersion law. This particular dispersion, that is only valid at 
low energies, mimics the physics of quantum electrodynamics (QED) for massless fermions except for 
the fact that in graphene the Dirac fermions move with a Fermi speed i>f, which is 300 times smaller 
than the speed of light c. Hence, many of the unusual properties of QED can show up in graphene but 
at much smaller speeds [HOE]. In addition Dirac fermions behave in unusual ways when compared 
to ordinary electrons if subjected to magnetic fields, leading to new physical phenomena [8] such as 
the anomalous quantum Hall effect, which was observed experimentally [21 [3]. 

On the other hand, many efforts have been employed to understand the scattering behavior of 
fermions in graphene, among them we cite for instance [9J [TOl 111) . In particular the authors in [9] 
calculated the mode-dependent transmission probability of massless Dirac fermions through an ideal 
strip of graphene (length I, width w, in absence of impurities or defects), to obtain the conductance 
and shot noise as a function of the Fermi energy. They found that the minimum conductance of order 
e 2 /h at the Dirac point (when the electron and hole excitations are degenerate) is associated with a 
maximum of the Fano factor (the ratio of noise power and mean current). For short and wide graphene 
strips the Fano factor at the Dirac point equals 1/3, three times smaller than for a Poisson process. 
This is the same value as for a disordered metal, which is counter intuitive since the dynamics of the 
Dirac fermions in graphene is ballistic. 

To generalize the analytical approach developed in [9] and study other type of scattering, we 
investigate the behavior of massless Dirac fermions in a flat chip of graphene. This will be based on 
the 2D massless Dirac equation which constitutes a good description of the low energy excitations of 
the original honeycomb lattice [12^ [T3l fH] . The flat chip is under the influence of a short range double 
barrier potential with infinite mass boundary condition [15], which then result in a quantization of 
the wave number associated with the confining y-direction. We consider the transmission of massless 
Dirac fermions through such a double barrier structure for various barrier parameters. Continuity of 
the wavefunction at each interface along with the infinite mass boundary condition in the y-direction 
resulted in a system of eight algebraic equations for eight unknown coefficients. The detailed study 
of such system gave rise to a variety of interesting situations that were investigated. In particular the 
single barrier results [5] were obtained as a particular case in our model. Then we studied the effect 
of different potential parameters on the transmission, conductance and shot noise. 

An interesting situation arises when we set the potential floor in the leads to zero, then our 2D 
problem reduces effectively to a ID massive Dirac equation with an effective mass proportional to 
the quantized wave number along the transverse direction called the y-direction. Thus confinement 
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along the y-direction generated effective masses for our fermions, which depend on the quantized wave 
number and its energy line spacing is proportional to the inverse of its width. On the other hand, we 
conclude that it is fairly interesting that the minimal conductivity and maximal Fano factor remain 
the same independently of the ratio between the two potentials barriers involved. 

The present paper is organized as follows, in section 2, we solve the Dirac equation to derive 
the energy spectrum using the infinite mass boundary condition along y-direction. These quantized 
wave numbers will serve as channel labels in dealing with tunneling in section 3 where the reflection 
and transmission amplitudes are determined through the electric current density. We then explore 
the obtained results in section 4 to discuss different situations related to the characteristics of the 
barriers. We have also devised a parameter which correlates the barrier heights and depth so as to 
enable us to investigate different potential configurations. The conductance and the Fano factor, for 
the present system, will be analyzed in section 5. In section 6, we show that the results obtained in [9] 
are a special case of our findings. We then conclude our work in the final section. 



2 Theoretical model 

In the system made of graphene, the two Fermi points, each with a two-fold band degeneracy, can 
be described by a low-energy continuum approximation with a four-component envelope wavefunction 
whose components are labeled by a Fermi-point pseudospin = ±1 and a sublattice forming a honey- 
comb. Specifically, the Hamiltonian for one-pseudospin component in the vicinity of the K point and 
in the presence of a scalar potential V{x) can be described by 

H = v F a ■p + V(x)I 2 (1) 

where the pseudospin matrices a are represented by the Pauli matrices, p = —ih V and Vf = 3ta/ (2ti) ~ 
10 6 ms _1 is the Fermi velocity of the massless Dirac fermions, t being the nearest neighbor hopping 
matrix element and a the carbon-carbon interatomic distance. The above equation can be written 
explicitly in matrix form using the system unit (h = c = e = Vf = 1) 

H = ( V(X) ~ idx " 9y ) (2) 

\ -id x + d y v(x) )' 

To go further, we fix the potential V(x) by considering a double barrier potential described picto- 
rially in Figure 1: 



d 2 - d 1 d 1 d 5 
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Figure 1: The double barrier potential along x-direction. 



This potential is defined by 

V , \x\ > d 2 

V(x) = I V u di < \x\ < d 2 (3) 

k V 2 , \x\ < d\ 

where V12 and d\ y2 are positive potential parameters (see Figure 1) characterizing the double barrier 
structure. We divide the configuration space into five regions numbered j = 1, • ■ • 5 associated with 
piece-wise constant potential sections. In the outside barrier regions V = Vq, V = V\ for j = 2, 4 in 
the barrier regions and finally V = V 2 for j = 3 in the well region. Each region j is characterized by its 
wave vector kj = kj x i + kjyj such that kj = k 2 x + fe| . The absence of potential along the y-direction 
leads to the conservation of kj y in all regions and its value is subsequently determined by the infinite 
mass boundary condition. 

Solving the eigenvalue equation 



V(x) -id x - d y 
-id x + d y V(x) 



i(j = Ei/j (4) 



we get plane wave solutions in all constant potential regions, such that in the j-th region we can write 

1 

V2 



^( X ,y) = J- I 1 )S-r (5) 



where the parameter Zj is given by 



^jx ^ r*jy 

with the sign function Sj = sign (E — Vj), the phase 6j = arctanf^J and the propagation wave 



number kj x = J (E — Vj) 2 — k 2 . The corresponding energy eigenvalue reads as 



E = Vj + Sjy/k 2 x + k] y (7) 

where the Sj = ±1 sign distinguishes electrons from hole conduction in the j-th constant potential 
region. One can show that the complex numbers Zj satisfies the identity 

1 + ^ = 2S 'P^% (8) 

which will play a crucial role in analyzing the tunneling effect and, in particular, in discussing different 
potential configurations used. 

One way to confine the particle along the transverse y-direction is to use a scalar potential (vector 
potential will not serve the purpose due the Klein tunneling) that rises to very large values at the 
boundary, the width of the strip along y-direction is denoted by w. This is achieved by using a 
mass term in the Dirac equation that rises to infinity at the edges y = and y = w. Otherwise, 
the infinite mass boundary condition gives us the quantification of the transverse wave vector kj y . 
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To determine this tangential component, we will consider a four-component eigenspinor, denoted by 
ij) = ifjn.kj( x ,y) = (V'l) V*2 5 V^)*- The first two components satisfy the Dirac equation without the 
mass term and the second two components satisfy the same equation by mapping pj y — > —pj y - For 
simplicity, we use the notations kj y = q n and kj x = kj to write the four-component eigenspinor as 

/ 1 \ / \ f zA I \ 



ip = a r , 





, , 



e i(q„y+kjx) _|_ 







e i(q n y+kjx) _|_ ^ 



V 1 / 



/, 3 \ 

1 


V J 



-i(q„y-kjx) , . / 




1 



-i(q n y-kjx) 



(9) 



The infinite mass boundary condition is expressed as follows [15 



a, 







y=0 



v -<7 X 

After some algebras, we end up with 



y=0, 



-a x 
a x 



y=W- 



7T / 1 

— \ n + - 

w\ 2 



(10) 



(11) 



where n is a positive integer number, which from now on denotes the propagation mode. Thereby 
confining the particle along the y-direction leads to a quantized value of the transverse wave vector. 
Finally, the energy spectrum is now given by 



(12) 



These results will be used to discuss different issues related to scattering phenomena, in particular, 
in the next section we will investigate the reflection and transmission coefficients through the double 
barrier structure. 



3 Transmission and reflection amplitudes 

The impact of the incident wave at a junction gives rise to reflected and transmitted waves. For the 
reflected wave, the wave vector kj along x-direction is opposite to that of the incident wave and the 
corresponding 9j angle is transformed into n — 6j. This allows us to write the spinors as 

^ n) {x, e *(M+^) (13) 




= ^[ _U ) e i( - kjX+qnV) (14) 



where the superscripts are obviously associated with waves traveling to the right (+) or to the left (— ) 
along the propagation direction (x). The subscripts, on the other hand, denote the scattering region 
(j) and the propagation mode (n). In the forthcoming analysis, we will see how these can be used to 
explicitly determine the reflection and transmission amplitudes. 

It is straightforward to solve the tunneling problem for Dirac fermions through the double barrier 
potential. Indeed, let us assume that the incident wave propagates at an angle 9\ with respect to the 
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positive x-direction and the impact of the incident wave on any junction gives rise to a reflected wave 
and transmitted wave. In this case, the upper and lower components of the Dirac spinor, denoted by 
ipi and 9321 axe given in each region by the following form 

( 



^gifcix _|_ re -ikix\ e iq n y 
( ae ik2X _|_ ] 3e -ik2X^ e iq n y 

(ce ik3X + de- ik3X )e iqnV ', 



ee 



ik2X 



+ fe 



-ik 2 x^ e iq n y 



^gifciXg iq„y 



^2 



(z±e lklX - 
(az2e zk2X 
(cz 3 e lk ' iX ■ 
(ez 2 e ik2X - 



rz^ l e~ iklX )e iqnV , 

bz^ l e- ik2x yi n y, 
dz^e-^y^y, 

f Z2 1 e ~ ik2 x ) e iq " y , 



t Zl e iklX e iqny , 



X < — G?2 
— C?2 < X < — d\ 

—d\ < x < d\ 
d\ < x < d 2 
x > d 2 

X < — C?2 

— C?2 < x < —d\ 
—d\ < x < d\ 
d\ < x < d 2 

X > G?2. 



(15) 



(16) 



Using the continuity of the spinor wavefunctions at the four potential discontinuities, one arrives at a 
system of eight algebraic equations for the eight unknown coefficients (a, b, c, d, e, /, r, t). These rela- 
tionships can be expressed in terms of transfer matrices between different regions but this sophisticated 
technique is not very much needed in our case. Solving this system of equations for the transmission 
amplitude, which is of interest to us, gives 

(17) 



B+C+D+E+F 



where the involved parameters read as 

A = eM[(te+*s)di-(*i-*a)«fe](l + ^)(l + ^)2(i + 



B 
C 
D 
E 
F 



Aik2d,2 



(Zl - z 2 ) 2 (z 2 - z 3 y 



3 4i ^ +k ^ dl (l + Z 1 Z 2 ) 2 (z2-Z 3 ) 2 
3 4i ^ dl+k2d2) ( Zl - Z 2 ) 2 (l + Z 2 Z 3 f 



(18) 



a 4ik 2 d 1 



(l + z 1 z 2 ) 2 (l + z 2 z 3 f 



2e 2ifc 2 (d 1 +d 2 ) ( _ 1 + g4*fc3di)( Zl _ Z2 ^ Z2 _ Z3 ) (1 + ZlZ2){l + 



Z2Z 3) 



To evaluate the reflection and transmission coefficients, we introduce the electric current density 
J for our system. After calculation, we obtain 



J = -iiiilvaij) 



(19) 



where ip stands for ip m = ifft n \(x, y), ^ ref = rtp^ n ^(x, y) and ip tr = tif)% n ){ x i V)- Computing explicitly 
(|19p gives for the incident, reflected and transmitted current density components 



Jx 



±i(zi + z*) = ±2is\ 



\Jk\ + 



J 



ref 



J 



tr 



=pir*r(zi + z*) = =p2ir*rsi 
±it*t(z 5 + z%) = ±2it*ts 5 



fci 



V k 't + il 
h 



(20) 



5 



The transmission and reflection coefficients, are expressed as follows 



141 

141 



Vg + g | t |2 



\K\ \t\ 



R 



J 



ref I 



(21) 
(22) 



In our case due to the symmetry of the potential configuration in the incident and transmission regions 
we have \K\ = 1, i.e. k\ = k§, and therefore we have T = \t\ 2 . The above results will be investigated 
numerically for different potential configurations to enable us to extract more conclusions regarding 
the basic features of our system. 



4 Limiting cases and discussions 



Similarly to the single barrier problem analyzed in [IB], let us investigate the energy spectrum struc- 
tures to understand further our system. Recall that, the spectrum of Dirac fermions in single-layer 
graphene is linear at low Fermi energies. Applying a potential barrier Vj in region j causes a displace- 
ment of the spectrum by an offset of Vj as shown in Figure 2. Each conical spectrum of graphene 
is the result of intersection of the energy bands originating from sublattices A (shown in red) and B 
(shown in blue). The portion of the spectrum, which is above the offset, represents the conduction 
band and the lower one represents the valence band. 






Figure 2: Tunneling through double barriers Vj in graphene where (a): Vq < E < V2, (b): V2 < E < V\, 
(c): E > V\. The green filling indicates occupied states and the pseudospin a is parallel (antiparallel) to 
the direction of motion of electrons (holes). 



The value of energy E compared to potential barrier strength Vj fixes the value of signature Sj = ±1 
in each region j. This sign plays an important role because it determines the nature of the particles 
(electrons or holes) involved in a given region j as well as the direction of wave vectors. Clearly, the 
symmetry between the conduction and valence bands allows us to take E > and the other case can 
be treated in similar way. 
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We are now in a position to consider solutions corresponding to different physical potential pa- 
rameters and find the associated transmission coefficients. First we start by treating the situation 
that generalizes [9] by considering V\ = V and Vi = a V, a being a constant parameter. The above 
potential parametrization will enable us to consider a variety of interesting physical situations and 
therefore simplify the analysis of the transmission (|17p . First let us consider the case Vo = and 
a = 1, which represents a single barrier potential problem. In this case (|17p becomes 

ln e 2 ld k 2{zi _ Z2) 2 + {1 + ZlZ2) 2 W 

where d = 2di = Ad\. Here, tunneling exists only for some allowed energies, which depend essentially 
on the quantized wave vector q n along y-direction. In Figure 3 we plot the transmission coefficient as 
a function of energy for different values of the potential parameters. 




(c): V < 0, with V = V 2 , V = 0, a = 1. 

Clearly, in Figure 3a, one can notice that there are four energy zones that characterize the transmission 
coefficient: 

• The first zone is bounded by the energy interval < E < q n , which can be seen as a forbidden 
zone. 
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• The second is the Klein zone, i.e. q n ^ E ^ V — <2n where V is the step/barrier height cine! q n the 
dynamical particle mass, which contains oscillations (resonances) in the transmission coefficient. 
This is the situation in which only oscillatory solutions exist throughout and where the so called 
Klein paradox reigns [18] . The physical essence of the Klein tunneling lies in the prediction that 
according to the Dirac equation, fermions can pass through strong repulsive potentials without 
the exponential damping expected in usual quantum tunneling processes. The energy range in 
which this oscillatory behavior persists is called the Klein zone 

• The third zone V — q n < E < V + q n is a bowl (window) of zero transmission. 

• Finally, the fourth zone E > V + q n contains the usual high energy barrier oscillations and 
asymptotically goes to unity at high energy. 

On the other hand, compared to our previous work |17j . we found a strong correlation between the 
present 2D massless Dirac fermions and massive ID Dirac fermions. That is, our 2D system with 
m = is equivalent to ID system with an effective mass ?n* — q n . Our effective ID system behaves as 
if its carriers have a dynamical mass that depends on the transverse quantized wave number. We can 
say that our effective ID system has many carriers each associated with a quantized effective mass 
depending on the transverse propagation mode n. Furthermore, according to Figure 3a and 3b one 
can draw interesting conclusions. Indeed, if m* < ^ then the transmission is zero for E < m* while 
for m* > the transmission vanishes E < V + m*. In addition, the isolated peak (T <C 1) in Figure 
3b is resulted from the fact that the allowed energy started from the value E = and specifically it 
appeared at the value E = 4 which corresponds to V = 8. However if the potential V is negative, 
we end up with a potential well behavior where the corresponding transmission is plotted in Figure 
3c. Note that for any pair of (m*,V) the allowed energy verifies the condition E > m*. In summary 
we can claim that the infinite mass boundary condition generated a dynamical mass for our original 
massless system, which is equivalent to space compactification of graphene from 2D to ID |17j . 

Let us treat the double barrier case where Vq = and a ^ 1, i.e. V% < V\ and V% > V\. In both 
cases, the transmission coefficient is plotted in Figure 4: 




Figure 4: Transmission coefficient as a function of energy for (a): {Vq = 0, a < 1) and (b): (Vo = 0, a > 1). 

In Figure 4a we distinguish six different zones characterizing the behavior of the transmission coeffi- 
cient. These are 



8 



• The first is a forbidden zone where < E < m*. 

• The second is the lower Klein energy zone characterized by resonances and m* < E < Vi — m* . 
Here we have full transmission at some specific energies despite the fact that the particle energy 
is less than the height of the barrier. Note that this happens only if V > 2m* . At these energies 
the wavefunction inside the barrier is oscillatory. 

• The third zone V2 — m* < E < V2 + m* is a window where the transmission is zero, the 
wavefunction is damped and transmission decays exponentially. 

• The fourth V2 + m* < E < V\ —m* is the higher Klein energy zone with transmission resonances. 

• The fifth V\ — m* < E < V\ + m* is a window where the transmission is mostly zero but contains 
resonance peaks corresponding to the bound states associated with the double barrier. Usually 
they are associated with the eigenvalues of the associated Hamiltonian or poles of the associated 
Green function. 

• The sixth zone E > V\ + m* contains oscillations, the transmission converges to unity at high 
energies similarly to the non-relativistic result. 

Noting that Figure 4a for massless 2D Dirac fermions coincides with transmission coefficient of massive 
2D Dirac fermions |19j . Contrary to the case a > 1 where V2 becomes larger than V\ and therefore 
the results are modified, see Figure 4b. Indeed, compared to Figure 4a, the behaviors of some zones 
are completely reversed like for instance the window zone. 

One can also consider the interesting case Vq = V2 = 0, which is experimentally easier to real- 
ize than Vq = —00 and V2 7^ 0. Figure 5 shows the behavior of the transmission coefficient, which is 
completely changed and shows only one Klein zone compared to Figure 4 which shows two Klein zones. 




Figure 5: Transmission coefficient as a function of energy for (a): (Vq = V2 = 0, V\ = 4) and (b): 

(V = V 2 = 0,V 1 = 8). 

At this stage, let us consider the case where Vq — > —00 and focus only on a > 1. In this situation, 
the transmission amplitude for the n-th mode reduces to 

2 e 2i[(k2+k 3 )d 1 -(k 1 -k 2 )d 2 ]^i + z\f(\ + z f) 

tn = A' + e 4i ( fc 3 rf i+ fc 2 rf 2)(l - z 2 ) 2 (l + z 2 z 3 ) 2 + e 4ifc 2 rf i(l + z 2 ) 2 {l + z 2 z 3 ) 2 ^ 
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where the parameter A' is given by 

A' = e 4ik ^(z 2 -l) 2 (z 2 -z 3 ) 2 + e^ k2+k ^(l + z 2 )\z 2 -z 3 ) 2 

+2e 2 i fc 2 (d 1+ d 2 ) ( _ 1 + e 4ifc 3 d 1)(1 _ z 2^ 2 _ Z3)(1 + ^ Z3 y (25) 

It is clear that the transmission coefficient depends on many parameters, such as the dimensions of the 
graphene sheet (w,l), the potential parameters (di,d 2 , V, a), the propagation mode n and the energy. 
To study this behavior we proceed by varying the potential parameters w, I and n while E and a are 
fixed. Figure 6 shows that the transmission coefficient contains two bowls, one centered at V = E 
with resonant peaks of unit amplitude and the other at V = — with vanishingly small transmission. 
The geometry of the graphene system has different effects on the transmission coefficient. Indeed, 
when w increases the bowl width decrease and at certain value of w the bowl depth starts to decrease. 
However, if I increases the number of oscillations of transmission increases as well. 




Figure 6: Illustration of the influence of each parameter of (n,w,l) on the transmission coefficient as a 
function of potential V. 

5 Conductance and Fano factor 

To complete our study, we investigate two important physical quantities relevant to double barrier 
structure in graphene, the conductance and the Fano factor. At zero temperature, they can be 
expressed, respectively, as 

G = 50 £ T - F = ^ n =°ti T } (26) 

n=0 l^n=0 ^ n 
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where T n = \t n \ 2 is the transmission coefficient for n-th mode, go = ^f- and the factor 4 accounts for 
the spin and valley degeneracy. Taking into account the fact that w 3> I, we plot the conductance 
a = G and Fano factor F versus the potential for different values of the energy. This resulted in a 
similar graph to that in [9] but shifted by an amount E along the potential axis. This is an obvious 
result since the at zero temperature E represent the Fermi level. 




V 



-10 -5 5 10 15 20 -5 5 10 15 



Figure 7: The conductance a and Fano factor F as a function of the potential V for different values of 
energy, with Vq = — oo,a = 1. 

One can see that the conductance and the Fano factor keep the same generic behavior for all energies 
except that they have been translated along potential axis, which means both of them are energy in- 
dependent. In the Figure 8 we show the behavior of the conductance and the Fano factor as function 
of the potential V different values of a. From this figure it clear that a possess a minimum value 
(7 = ^ but F has a maximum value F = | at zero energy independently of the value of a. This result 
generalizes that obtained in [9], which was limited to the special case a = 1 




Figure 8: The conductance a and Fano factor F as a function of V at the Dirac point E = 0, for dif- 
ferent values of a (blue curve for a = 1, red curve for < a < 1 and green curve for a > 1), with Vq = — oo. 



11 



From Figure 8 we can conclude that for all values of a, we have the same conductance minimum a = ^ 
and Fano factor maximum F = ^ at the Dirac point (E = 0) but two extra satellite minima/maxima 
have been born. 

It is worthwhile to investigate what happens if the energy is E ^ and a > 1. This is described in 
Figure 9, which shows two conductance minima at the points E and — and similarly for the two max- 
ima that occur for the Fano factor. It is interesting to note that the first minimum in the conductance 
(Figure 9a) has no internal structure while the second one seems to have an oscillatory modulation 
reminiscent of some interference phenomena. A similar observation can be made on Figure 9b for the 
Fano factor. 




Figure 9: The conductance a and Fano factor F in terms of the potential V for E > and a = 4, with 

Vo = — oo. 

It is interesting to underline the behavior of the conductance a and Fano factor F in terms of the 
inter-barrier distance d\ for V2 = 0. The present case is plotted for different values of the potential V 
in Figure 10: 
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Figure 10: The conductance a and Fano factor F versus the inter-barrier distance d\ for V 2 = 0, with 

V = -00. 



As shown in Figure 10, G is decreasing as long as d\ is increasing and it is going to hold a constant 
value for large d\. However this is not the case for F, which is increasing with d\ up to certain value 
and becomes almost constant for a large interval of d±. As it is clearly shown in Figure 10b and lOd, 
the behavior of G and F are not affected so much by the variation of the potential V > 1. 



6 Recovering the single barrier results 

Let us show how to recover one single barrier results presented in [9] from what we obtained and 
discussed so far. We start by recalling that in [S] the authors considered the single barrier case with 
the potentials Vq — > — 00 and V\ = Vi = V. This allows us to end up with the following restrictions 
on our parameters 

si = s 5 = 1, ki = k 5 = 00, Zl = z 5 = l (27) 
S2 = s 3 = s 4 , k 2 = h = h, z 2 = z 3 = z 4 . (28) 



Injecting the above results in (|17p and f)18[) . we get the transmission amplitude 

2 e - 2 ^(fci-te)( 1 + 2) 
tn ~ e 4id ^(l-z 2 ) 2 + (l + z 2 ) 2 ' 



(29) 



Since the potential Vo is negative then the "wave vector k\ is purely recti ctnd. consequently t n can be 
reduced to the form 

2(1 + z$) 

tn ~ e ilk 2(l-z 2 ) 2 + e- ilk 2(l + z 2 ) 2 (30) 
where I = 2d 2 = 4d±. Figure 11 shows the transmission behavior versus V for two values of propaga- 
tion mode. 
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Figure 11: Transmissions coefficient as a function of the potential V for two values of the propagation 
mode n = 4 and n = 5, with Vq = — oo. 



It is clear that two transmission curves are symmetric with respect to the point V = E while an 
increase in the quantum number n widen the bowl width. On the other hand, setting Vo —> — oo and 
V — > 0, the transmission probability T n associated with the n-th mode at the Fermi level is obtained 
from T n = \t n \ 2 = \t\ 2 for E = 0. Under these circumstances our parameters are given by 



h = iq n , z 2 
then our transmission amplitude takes a simple form 



oo 



2e lq " 



\ _|_ e 2lq„ 



cosh(lq n 



cosh 



~irl 






H)] 


w 





(31) 



(32) 



which coincides with that of [9]. Using (|26p . it is shown that a minimum in the conductance and a 
maximum in the Fano factor [9] occur at the Dirac point. In the limiting case, i.e. j — > oo (for a 
short and wide strip) at the Dirac point a and F, respectively, reduce to 



9o j-, 1 
(7 ~ — , r ~ — . 

7T 3 



(33) 



Now let us consider a nonzero potential V and evaluate the transmission coefficient at the Dirac 
point E = to see how the transmission obtained in [9] will behave. Indeed, from our previous result 
dSJ) we can easily show that 



1 + zi 



lTVk 2 . 



Using these results we obtain transmission coefficient T n for the n-th mode as follows 

2 



(12 



k,2 cos(k2l) + iV sin(£>20 



(34) 



(35) 



To illustrate this case and identify the difference with respect to Figure 11, we plot the above result 
for two particular values of n in Figure 12: 
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Figure 12: Transmissions coefficient as a function of V at the Dirac point (E = 0) for n = 4 and 
n = 5, with Vo = — oo. 

One can see that both curves start from zero transmission and oscillate to reach a total transmission, 
the valley gets wider as n increases. This behavior is similar to the ID massive Dirac equation if 
one increase the mass of the Dirac Fermion, hence establishing a strong relationship between the two 
systems. 

7 Conclusion 

We have generalized the results obtained in [9] to double barrier case. More precisely, we have analyzed 
transport properties through double barrier structure in graphene. In the beginning, we have solved 
the Dirac equation, for different regions forming the whole system made of graphene, to get the 
solutions of the energy spectrum. In the present analysis, the infinite mass boundary condition was 
taken into account, which then resulted in a quantization of the wave number associated with the 
confining y-direction. This later is interpreted as a dynamical effective mass of Dirac fermions and 
compared to previously works. 

Subsequently, the transmission of massless Dirac fermion through the double barrier structure for 
various barrier parameters is studied. Continuity of the wavefunction at each interface along with the 
infinite mass boundary condition in the y-direction resulted in a system of eight algebraic equations 
for eight unknown coefficients. The detailed study of such system gave rise to a variety of interesting 
situations that were investigated. 

In fact, we have studied the effect of different potential parameters on the transmission, conduc- 
tance and shot noise. An interesting situation arises when we have set the potential floor in the leads 
to zero, then our 2D problem reduces effectively to a ID massive Dirac equation with an effective 
mass proportional to the quantized wave number along the transverse y-direction. Thus confinement 
along the y-direction generated effective masses for our fermions, which depend on the quantized wave 
number and its energy line spacing is proportional to the inverse of its width. Furthermore, we noticed 
that the minimal conductivity and maximal Fano factor remain the same independently to the ratio 
between the two potentials ( = ol ) • 



Finally, we have shown how to recover from our model the single barrier results [9j. This was 
done by considering the following limit Vq = and V\ = V2 = V in our general formulation. For the 




15 



single barrier case with Vo = — oo and V\ = V2 = V we have recovered the result of reference [9]. This 
showed that our results are interesting an more general. 
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